
Some exact asymptotics 
in the counting of walks 
in the quarter-plane 

Guy Fayolle, Kilian Raschel 







z 




RESEARCH 


z 

DC 
CO 




REPORT 


cn 




N° 7863 


CO 

to 




January 2012 


CM 

o 
Z 




Project-Team Imara 


co 
co 




Some exact asymptotics in the counting of 
walks in the quarter-plane 

Guy Fayolle*, Kilian Raschel t 
Project- Team Imara 
Research Report n° 7863 — January 2012 — 17 pages 



Abstract: Enumeration of planar lattice walks is a classical topic in combinatorics, at 
the cross-roads of several domains (e.g., probability, statistical physics, computer science). 
The aim of this paper is to propose a new approach to obtain some exact asymptotics for 
walks confined to the quarter plane. 
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Quelques estimations asymptotiques exactes pour le 
comptage de marches dans le quart de plan 



Resume : L 'enumeration de marches planaires sur des treillis est un domaine classique 
de la combinatoire, a la croisee de plusieurs disiplines (probabilites, physique statistique, 
informatique) . Cet article propose une nouvelle methode donnant des estimations asymp- 
totiques exactes pour certaines marches confinees dans le quart de plan. 

Mots-cles : Analyse de singularity, fonction generatrice, marche aleatoire dans le quart 
de plan, probleme aux limites. 
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Fig. 1.1: Four famous examples, known as the simple, Kreweras', Gessel's and Gouyou- 
Beauchamps' walks, respectively 



1 Introduction 

Enumeration of planar lattice walks is a most classical topic in combinatorics. For a given 
set S of admissible steps (or jumps), it is a matter of counting the number of paths of a 
certain length, which start and end at some arbitrary points, and might even be restricted 
to some region of the plane. Then three natural important questions arise. 

(i) How many such paths do exist? 

(ii) What is the asymptotic behavior, as their length goes to infinity, of the number of 
walks ending at some given point or domain (for instance one axis)? 

(iii) What is the nature of the generating function of the numbers of walks? Is it holo- 
nomic§, and, in that case, algebraic or even rational? 

If the paths are not restricted to a region, or if they are constrained to remain in a half- 
plane, it turns out [6] that the counting function has an explicit form (question (i)), and 
is, respectively, rational or algebraic (question (iii)). Question (ii) can then be solved from 
the answer to (i). 

The situation happens to be much richer if the walks are confined to the quarter plane 
Z+. As an illustration, let us recall that some walks admit an algebraic counting function, 
see [14] for Kreweras' walk (see Figure 1.1), while others admit a counting function which 
is not even holonomic, see [6] for the so-called knight walk. 

In the sequel, we focus on walks confined to Z+, starting at the origin and having small 
steps. This means exactly that the set S of admissible steps is included in the set of the 
eight nearest neig hbors, i.e., S C {-1, 0, l} 2 \ {(0, 0)}. By using the classical Kronecker's 
delta, we shall write 

1J \ if $ s. y ' 

On the boundary of the quarter plane, allowed jumps are the natural ones: steps that 
would take the walk out Z+ are obviously discarded, see examples on Figure 1.1. 
There are 2 8 such models. Starting from the existence of simple geometrical symmetries, 
Bousquet-Melou and Mishna [5] have shown that there are in fact 79 types of essentially 
distinct walks — we will often refer to these 79 walks tabulated in [5]. For each of these 
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79 models, q(i,j,n) will denote the number of paths confined to 71?+, starting at (0,0) and 
ending at (i,j) after n steps. The associated generating function will be written 

Q{x,y,z) = q(i,j,n)x i y j z n . (1.2) 

i ,j,n^0 

Answers to questions (i) and (iii) have been recently determined for all 79 models of walks. 
To present these results, we need to define a certain group, introduced in a probabilistic 
context in [18], and called the group of the walk. This group W = of birational 

transformations in C 2 leaves invariant the function ^-i<i j<i ^i,j x% V^ -> an< ^ nas the- two 
following generators: 

*<*,,) = te-ssi^iA (i.3) 

Clearly x I'o^/ = $o$ = id, and W is a dihedral group, whose order is even and at least 

four. The order of W is calculated in [5] for each of the 79 cases: 23 walks admit a finite 

group (then of order 4,6 or 8), and the 56 others have an infinite group. 

The expression of the counting function Q(x, y, z) — or equivalently the answer to (i) — for 

the 23 walks with a finite group has been determined in [2, 5], It was basically derived 

by applying a kernel method (see, e.g., [13]) to the functional equation (2.1) satisfied by 

Q(x,y,z). 

As for the 56 models with an infinite group, the function Q(x, y, z) related to the 5 singular 
walks (i.e., walks having no jumps to the West, South- West and South) was found in 
[19, 20] by using elementary manipulations on (2.1). For the remaining 51 non-singular 
walks, Q(x,y,z) was found in [21] by the method, initiated in [8, 9], of reduction to 
boundary value problems (BVP). 

Concerning (iii), the answer was obtained in [2, 5, 10] for the 23 walks with a finite group: 
the function Q(x,y,z) is always holonomic, and even algebraic when ijSij is 

positive. For the 56 models with an infinite group, it was also proved that they all admit 
a non-holonomic counting function: see [19, 20] for the 5 singular walks, and [16] for the 
non-singular ones. 

Let us now focus on question (ii), which actually is the main topic of this paper. A 
priori, it has a link with question (i) (indeed, being provided with an expression for the 
generating function, one should hopefully be able to deduce asymptotics of its coefficients), 
and with (iii) as well (since holonomy is strongly related to singularities, and hence to the 
asymptotics of the coefficients). As we shall see, the situation is not so simple, notably 
because the expression of the counting function (1.2) is fairly complicated to deal with, in 
particular when it is given under an integral form, which is a frequent situation in concrete 
case studies. 



§ A function of several complex variables is said to be holonomic if the vector space over the field of 
rational functions spanned by the set of all derivatives is finite dimensional. In the case of one variable, this 
is tantamount to saying that the function is solution of a linear differential equation where the coefficients 
are rational functions (see [13]). 
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If we restrict question (ii) to the computation of asymptotic estimates for the number of 
walks ending at the origin (0, 0) (these numbers are the coefficients of Q(0, 0, z), see (1.2)), 
a very recent paper [7] solved the problem for much more general random walks in M. d , for 
any d > 1: indeed, the authors obtain the leading term of the estimate, up to a coefficient 
involving two unknown harmonic functions. On the other hand, it is impossible to deduce 
from [7] the asymptotic behavior of the number of walks ending, for instance, at one axis 
(coefficients of Q(1,0, z) and Q(0,l,z)), or anywhere (coefficients of Q(l,l,z)). 
Still with respect to question (ii), it is worth mentioning several interesting conjectures by 
Bostan and Kauers [3, 4], which are proposed independently of the finiteness of the group. 
The aim of this paper is to propose a uniform approach to answer question (ii), which, as 
we shall see, works for both finite or infinite groups, and for walks not necessarily restricted 
to have closed trajectories as in [7]. 

First, Section 2 presents the basic functional equation, which is the keystone of the analysis. 
Then, Section 3 illustrates our approach by considering the simple walk (the jumps of 
which are represented in Figure 1.1), which has a group of order 4. In this case, we obtain 
the explicit expression of Q(x,y,z), allowing to get exact asymptotics of the coefficients 
of Q(0, 0, z), Q(l, 0, z) (and by an obvious symmetry, also of Q(0, 1, z)) and of Q(l, 1, z), 
by making a direct singularity analysis. Section 4 sketches the main differences (in the 
analytic treatment) occurring between the simple walk and any of the 78 other models. In 
particular, we explain the key phenomena leading to singularities. Indeed, a singularity is 
due to the fact that a certain Riemann surface sees its genus passing from 1 to 0, while 
another singularity takes place at z = 1/|<S|, where |<S| denotes the number of possible 
steps of the model. Also, as in the probabilistic context (although this will not be shown 
in this paper), the orientation of the drift vector has a clear impact on the nature of the 
singularities which lead to the asymptotics of interest. 

To conclude (!) this introduction, let us mention that our work has clearly direct con- 
nections with probability theory. Indeed, with regard to random walks evolving in the 
quarter plane, and more generally in cones, the question of the asymptotic tail distribu- 
tion of the first hitting time of the boundary has been for many decades of great interest. 
The approach proposed in this paper should hopefully lead to a complete solution of this 
problem for random walks with jumps to the eight nearest neighbors. 

2 The basic functional equation 

A common starting point to study the 79 walks presented in the introduction is to estab- 
lish the following functional equation (proved in [5]), satisfied by the generating function 
defined in (1.2). 

K(x,y,z)Q(x,y,z) = c(x)Q(x,0,z) +c(y)Q(0,y,z) - £_i ,-iQ(0, 0, z) - xy/z, (2.1) 
where c and c are defined in (2.3) and 

K(x,y,z) = xyE-Kij^i ki xi V j ~ V*]- ( 2 - 2 ) 
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Equation (2.1) holds at least in {|a;| ^ 1, \y\ ^ 1, \z\ < 1/|<S|}, since obviously q(i,j,n) ^ 
|5| n . In (2.1), we note that the <$ij's defined in (1.1) play the role (up to a normalizing 
condition) of the usual transition probabilities Pi j's in a probabilistic context. 
Our main goal is to analyze the dependency of Q(x, y, z) with respect to the time variable 
z, which in fact plays merely the role of a parameter, as far as the functional equation 
(2.1) is concerned. Remark that Q(x,y,z) has real positive coefficients in its power series 
expansion. Thus, letting R denote the radius of convergence, say of Q(1,0, z), it is worth 
noting that, by Pringsheim's theorem (see [23]), z = R is a singular point of this function. 
The quantity (2.2) that appears in (2.1) is called the kernel of the walk. It can be rewritten 
as 

a(y)x 2 + [b(y) - y/z]x + c(y) = a(x)y 2 + [b(x) - x/z]y + c(x), 

where 



a ( x ) = E-ia<i s i,i xi > = E-i<i<i Siflx*, c(x) = E-i<i<i k-ix*. 



(2.3) 



Let us also define 

d(y, z) = \b{y) - y/z} 2 - Aa{y)c{y), d{x, z) = [b(x) - x/z] 2 - Aa{x)c{x). (2.4) 

From now on, we shall take z to be a real variable. Indeed this is in noway a restriction, 
as it will emerge from an analytic continuation argument. 

For z G (0, 1/|<S|), the polynomial d (resp. d) has four roots (one at most being possibly 
infinite) satisfying in the x-plane (resp. y-plane) 

\xi(z)\ < x 2 (z) < 1 < x 3 (z) < |x 4 (z)| < oo, \yi(z)\ < y 2 (z) < 1 < y 3 (z) < \y&(z)\ «S oo, 

(2.5) 

as shown in [21]. Then consider the algebraic functions X(y,z) and Y(x,z) defined by 

K(X(y, z),y, z) = K(x, Y(x, z), z) = 0. 
With the notations (2.3) and (2.4), we have 



)_-*(»■»> fV^foj, Ylx , z) = zM^lp/^A. (2 . 6) 

2a(y,z) 2a(x,z) 

From now on, we shall denote by Xq, X\ (resp. Yq, Y\) the two branches of these algebraic 
functions. They can be separated (see [9]), to ensure \Xq\ ^ |^i|, resp. |Yo| ^ |li|, in the 
whole complex plane C. 

3 The example of the simple walk 

To illustrate the method announced in the introduction, we consider in this section the 
simple walk (i.e., 5ij = 1 only for couples such that ij = 0, see Figure 1.1). Then 
the main following results hold. 
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Proposition 3.1. For the simple walk, 



1 f 1 1 - 2uz - J (I - 2uz) 2 - 4z 2 I r , x 

gO,0,z =" Vl-u 2 du. 3.1 

7T J_ x Z l 

Note that Q(0, 0, z) counts the numbers of excursions starting from (0,0) and returning 
to (0,0). 

Proposition 3.2. For the simple walk, as n — > oo, 



g(0,0,2n) 
Proposition 3.3. For the simple walk, 



4 16" 
7r n 3 



1 f l \-2uz- J(\ - 2uz) 2 - 4z 2 ll + u , 

The generating function Q(l, 0, z) depicts the numbers of walks starting from (0,0) and 
ending when they reach the horizontal axis. In addition, by an evident symmetry, Q(0, 1, z) - 
Q{l,0,z). 

Proposition 3.4. For the simple walk, as n — > oo, 

/ \ 8 4" 

Having the expressions of Q(l, 0, z) and Q(0, 1, z), the series Q(l, 1, z) is directly obtained 
from (2.1). The asymptotics of its coefficients, which represent the total number of walks 
of a given length, is the subject of the next result. 
Proposition 3.5. For the simple walk, as n — > oo, 

V q{i,3,n) ~ . 

t—' it n 

i,j>0 

Proof of Proposition 3.1. The key point is to reduce the computation of Q(x, 0, z) to a 
BVP set on T = {t G C : \t\ = 1}, according to the procedure proposed in [8, 9]. Indeed, 
for any t £ T, we have 

c(t)Q(t, 0, z) - c(t)Q(t, 0, z) = tY o(t,z)-tY (t,z) ^ 

z 

where t stands for the complex conjugate of t. The proof of (3.2) relies on simple ma- 
nipulations on the functional equation (2.1). The main argument is the following. On 
K(x,y,z) = 0, letting y tend successively from above (y + ) and from below (y~) to an 
arbitrary point y on the cut [yi(z), 3/2(2)] defined in (2.5), Q(0,y,z) remains continuous, 
and thus can be eliminated by computing the difference Q(0,y + , z) — Q(0,y~ , z) = in 
the functional equation (2.1), see [8, 9] for full details. 
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In the present case, a pleasant point (from a computational point of view) is that, on the 
unit circle, we have t = 1/t. Hence, after mutiplying both sides of (3.2) by l/(t — x) with 
\x\ < 1, integrating over T, and making use of Cauchy's formula, we get the following 
integral form (noting that here c(x) = x) 

c(x)Q(x, 0, z) = J- f ^MzM!) d i, VM < 1. (3.3) 
2%-kz Jy t — x 

For t E T, we know from general results [9, 21] that for z E [0, 1/4], Yo(t,z) is real and 
belongs to the segment [yi(z),y2(z)], the extremities of which are the two branch points 
located inside the unit circle in the y-plane, see Section 2. Then it is not difficult to check 
that the integral in the right-hand side of (3.3) does vanish at x = 0. Hence we can 
entitled to write 

Q(0, 0, *) = hm -J— / tY (t,z)-jy (i,z) dt = J_f tY (t,z)- tY (i,z) 

^ v ' ' ; x^lmxz J r t-x 2mz J r t 2 ' v ' 

where the last inequality follows from l'Hospital's rule. 

To work on formula (3.4), it is convenient to make the straightforward change of variable 
t = e l9 , which yields 

Q(0,0,z) = — f sm9e- i9 Y (e i9 ,z)d6 = — [ sin 2 9 Y {e ie , z) dO, (3.5) 

2-KZ Jq 7TZ J 

where we have used the fact that Yo(e te , z) is real and Yo(e l9 , z) = Yo(e~ ie , z). Furthermore, 
we have by (2.6) 



Me", ,) = L^g^ySZjj^gEjg. (3.6) 
Then, instantiating (3.6) in (3.5), we get after some light algebra 

Q^Q jZ ) = J--J-f sin 2 6y/(l - 2z cos 6) 2 - Az 2 6.6. (3.7) 

2z TTZ Jq 

Putting now u = cosO in the integrand, equation (3.7) becomes exactly (3.1). ■ 

Proof of Proposition 3.2. The function Q{0, 0, z) is holomorphic in C \ ((— oo, — 1/4] U 
[1/4, oo)), as it easily emerges from Proposition 3.1. Accordingly, the asymptotics of 
its coefficients will be derived from the behavior of the function in the neighborhood of 
±1/4 (see, e.g., [13]). On the other hand Q(0, 0, z) is even, as can be been in (3.1), or 
directly since (/(0, 0, 2n + 1) = 0, for any n ^ 0, in the case of the simple walk (see Figure 
1.1). Consequently, it suffices to focus on the point 1/4. 



First, we rewrite y^l - u 2 )[(l - 2uz) 2 - 4z 2 } as ^(1 - u)(l - 2(u + l)z)y/(l + u)(l - 2(u 
where the second radical admits the expansion £V fJ>ij(u — l) l (z — 1/4)- 7 . Therefore, 
for z in a neighborhood of 1/4, we have 



1 Z" 1 

Q(0,0,*) = , Vftj(z-1/4) J (u - 1)V(1 - - 2{u + l)z) du. (3.8) 

7T Z — / i 
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For any % 0, there exist two functions, say fa and gi, which are analytic at z = 1/4 and 
such that, for z in a neighborhood of 1/4, 

f {u - l)V(l-«)(l-2(« + l)z)d« = - l/4) 4+2 ln(l - 4z) + (3.9) 

Now, by inserting identity (3.9) (which will be proved right after the proof of Proposition 
3.2) into (3.8), we obtain the existence of a function, say g, analytic at z = 1/4 and such 
that 

16 

Q(0,0,z) = Mo,o(^ - l/4) 2 ln(l - 4z)[/ (l/4) + 0(z - 1/4)] + «/(*). (3.10) 

7T 

Furthermore, it is easy to prove /o(l/4) = 4\/2 and ^o,o = v2. With these values, a 
classical singularity analysis (see, e.g., [13, 23]) shows that, as n — > oo, the nth coefficient 
of the function in the right-hand side of (3.10) is equivalent to (16/7r)4 n /n 3 . It is then 
immediate to infer that 

. 32 4 2n 
g 0,0,2n 3, 

remarking in the latter quantity the factor 32 (and not 16 !), due to the parity of Q(0, 0, z) 
mentioned earlier. ■ 

Proof of (3.9). Let Aj(z) denote the integral in the left-hand side of (3.9). The change of 
variable u = 1/(42) — [(1 — 4z)/(4z)]v gives 

fl-AzY /• 1 + 8 */(l-4*) . , 

^(z) = [—-^\ V2i J {l-v)Wv 2 -Idv. 

Then, by letting v = cosht, we conclude that 

/1-4*V i— /■ cosh_1 ( 1 +8 2 /( 1 -^)) 
Ai(z) = ( ^ \ V2z J (l-cosht) 4 (-l + cosh 2 t)dt. (3.11) 

Besides, by a classical linearization argument, there exist real coefficients cto, . . . , such 
that 

i+2 

(1 - cosht)*(-l + cosh 2 t) = Y^ «fc cosh(kt). (3.12) 

So, by means of (3.12), it is easily to integrate Ai(z). Then a delinearization argument 
shows the existence of two functions, gi and hi, which are analytic at z = 1/4 and satisfy 

Ai(z) = 9i (z) + 0-^-) v / 2i/ li (z)cosh- 1 (l + 8z/(l-4z)). 

Finally, remembering that, for all u > 1, cosh -1 u = log (it + Vu 2 — 1), equation (3.9) 
follows, since cosh -1 (l + 8z/(l — 4z)) + ln(l/4 — z) is analytic at z = 1/4. ■ 
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Proof of Proposition 3.3. The argument mimics the one used in Proposition 3.2. Instan- 
tiating x = 1 in equation (3.3) and using (3.6) yields directly 

= - j / (1 + cos 0)(1 — 2z cos 9 - v / (l-2zcos6») 2 -4z 2 ) 66, 

2-kz z J 

and the proof of Proposition 3.3 is terminated. ■ 

Proof of Proposition 3.4- It is quite similar to that of Proposition 3.2, by starting from 
the integral formulation of Q(l, 0, z) written in Proposition 3.3, and making an expansion 
of the integrand in the neighborhood of u = 1 and z = 1/4. So we omit it. ■ 

Proof of Proposition 3.5. It is an immediate consequence of Proposition 3.4, by using the 
equality 

(4-l/z)Q(l,l,z)=2Q(l,0,z)-l/z, 
which follows from (2.1). ■ 



4 Asymptotics of the number of walks with small steps confined to the 
quarter plane: a general approch 

This part aims at extending the approach of Section 3 to all 79 models. For this, we show 
in Section 4.1 that Q(x, 0, z) and Q(0, y, z) satisfy some BVP, that we solve in Section 4.2. 
Next, in Section 4.3 we compute Q(0, 0, z), Q(1,0, z), (5(0,1, z) and Q(l, l,z) (question 
(i)). In Section 4.4, we see that the analysis of Section 3 applies verbatim to the 19 
walks having a group of order 4. Finally, in Section 4.5 we analyze the singularities of the 
counting functions (question (ii)). In this short paper, we cannot provide the full details 
for the analysis of these singularities: we just explain where the singularities come from, 
and we postpone to the ongoing work [12] the computation of the exact behavior of the 
generating functions in the neighborhood of the (dominant) singularities. 



4.1 Reduction to a boundary value problem 

Most of the results in this section are borrowed from [9, 21]. In the general framework 
(i.e., for all 79 models), equation (3.2) holds on the curve, drawn in the x-complex plane 

C, ' 

M» = Xo{ [y 1 {z),y 2 {z% z) = X 1 ( [y 1 (z),y 2 (z)] , z), 

depicted in the next lemma, when the genus of the Riemann surface (corresponding to the 
manifold {(x,y :G C 2 : K(x,y,z) = 0}) is equal to 1. Here, remembering that z is real, 
we let 

[yi(z),y 2 (z)] 
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stand for a contour, which is the slit [yi(z), 7/2 (z)] traversed from yi(z) to 2/2(2) along the 
upper edge, and then back to yi(z) along the lower edge. 

Lemma 4.1. The curve M z is one of the two components of a plane quartic curve with 
the following properties. 

• It is symmetrical with respect to the real axis. 

• It is connected and closed in CD {00}. 

• It splits the plane into two connected domains, and we shall denote by ( £{M Z ) that 
containing the point x\(z). In addition, &(A4 Z ) C C \ [2:3(2:), £4(2)] . 

Similarly, we can define, in the y-complex plane C, the curve C z and the domain &(C Z ). 
However, in the case of a group (1.3) of order 4, the two components introduced above 
coincide to form a double circle, as in Section 3. Also, in the case of genus 0, there is only 
one component, which for instance can be an ellipse, see [9, Chapter 6]. 
Then, starting from the formal boundary condition (3.2) on the curve A4 Z , the function 
Q(x, 0, z) can be analytically continued with respect to x from the unit disc to the domain 
@(M Z ), see [9, 21]. The properties quoted in this section lead to the formulation of the 
fundamental BVP. 

Find a function x 1— > Q(x,0,z) analytic in the domain &(M Z ), and satisfying condition 
(3.2) on the boundary M. z . 

4.2 Solution of the boundary value problem by means of conformal gluing 

To solve the BVP stated in the previous section, we make use of [8, 9, 21]. In particular, 
we need the notion of conformal gluing function (CGF). Indeed, as it is, the BVP holds on 
a curve that splits the plane into two connected components (see Lemma 4.1). However, 
the BPV on the curve Ai z can be reduced to a BVP set on a segment, which is somehow 
computationally more efficient (see [17] and references therein). 

Definition 4.2. Let ?cCU {00} be an open and simply connected domain, symmetrical 
with respect to the real axis, and not equal to 0, C and C U {00}. A function w is said 
to be a conformal gluing function (CGF) for the domain & if it satisfies the following 
conditions. 

• w is meromorphic in c € . 

• w establishes a conformal mapping of ^ onto the complex plane cut along a segment. 

• For all t in the boundary of w(t) = w(J). 

For instance, the mapping 1 1— > t + l/t is a CGF for the unit disc centered at 0. It is worth 
noting that the existence (however without any explicit expression) of a CGF for a generic 
domain is ensured by general results on conformal gluing [17, Chapter 2]. In the sequel, 
we shall assume that the unique pole of w is at t = 0. 
The main result of this section is the following. 
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Proposition 4.3. For x G &(M Z ), 

c(x)Q(x,0,z)-c(0)Q(0,0,z) = -L / tY (t,z) W '^ Z \ -dt. (4.1) 

2mzJ Mz w(t,z) - w(x,z) 

Of course, a similar expression could be written for Q(0, y, z). 

4.3 Computation of Q(0, 0, 2), Q(l, 0, 2), Q(0, 1, z), Q(l, 1, 2) 

The expression we shall obtain for Q(0, 0, z) depends on c{x) or symmetrically on c(y), 
see equation (2.3), in the following respect. 

Suppose first c(0) = (this equality holds for the simple walk). Then, as in Section 3, we 
can write / 

Q(0,0,z) = lim f tY (t,z) . w ^ z \ dt. (4.2) 

x-*o 2mzc{x) J Mz w{t, z) — w[x, z) 

If c(0) 7^ and c{x) is not constant, then it has one or two roots, which necessarily are 
located on the unit circle V, see (2.3) and (1.1). Let x be one of these roots. Provided 
that x G &(Aig), we can write 

Q(0, 0, z) = / tY (t, z) W ' {t ' z) dt. (4.3) 

2mz J Mz w(t,z) - w(x,z) 

If c(0) 7^ and c(x) is constant, then we can evaluate the functional equation (2.1) at 
any point (x,y) such that \x\ ^ 1, \y\ ^ 1, and K(x,y,z) = to obtain an expression for 
Q(0,0,z). 

The computation of (3(1,0, z) and Q(0, l,z) depends on the position of the point 1 with 
respect to <£{M. Z ) and %?(£ z ). Indeed, if 1 belongs to these domains, then Q(1,0, z) and 
Q(0, l,z) are simply obtained by evaluating the integral formulations (4.1) at x = 1 and 
y = l. 

We now assume that for a given z, the point 1 does not belong to < £(A4 Z ). Then, by 
evaluating the functional equation (2.1) at (x,Yq(x, z)) and (Xq(Yq(x, z), z), Yq(x, z)), and 
by making the difference of the two resulting relations, we obtain 

c(x)Q(x, 0, z) = c(X (Y (x, z), z))Q(X (Y (x, z),z),0, z) + Y °^ Z \ X _ X (Y (x, z), z)]. 

(4.4) 

The key point is that, for x G &(M Z ), the range of C through the composite function 
Xq(Yo(x, z), z) is &(Ai z ) itself. This fact was proved in [9, Corollary 5.3.5] for z = 1/|<S|, 
but the line of argument easily extends to other values of z. In particular, to compute the 
right-hand side of (4.4), we can use the expression (4.1) valid for any x, in particular for 
x = 1. 

As for Q(l, 1, z), we simply use the functional equation (2.1), so that 

(|5| - l/z)Q(l, 1, z) = c(l)Q(l, 0, z) + c(l)Q(0, 1, z) - <5_i ,_iQ(0, 0, z) - l/z. (4.5) 
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4.4 Case of the group of order 4 

For the 19 walks with a group (1.3) of order 4 (the step sets S for these models have been 
classified in [5]), there are two possible ways to compute the exact expression and the 
asymptotic of the coefficients of Q(0, 0, z), Q(l, 0, z), Q(0, 1, z) and Q(l, 1, z). 
Firstly, we can use Section 4.2, as for any of the 79 models. Secondly, the reasoning 
presented in Section 3 extends immediately to the 19 walks having a group of order 4. 
Indeed, since in this case the boundary condition (3.2) is set on a circle, we can replace t 
by a simple linear fractional transform of t — for example, for the unit circle centered at 0, 
t = 1/t. In addition, it is proved in [9] that having a BVP set on a circle is equivalent for 
the group to be of order 4. 

For the simple walk, one can check that these two ways of approach coincide: take t + l/t 
for the CGF w(t,z) in (4.1), then make a partial fraction expansion of w'(t,z)/[w(t,z) — 
w(x,z)], and this will eventually lead to (3.3). 

4.5 Singularities of the generating functions 

As written in the introduction of Section 4, we cannot in this short paper go into deeper 
detail about the analysis of singularities. Let us simply state that only the real singularities 
of Q(0, 0, z), Q(l, 0, z), Q(0,1, z) and Q(l,l,z) play a role. Below, from the expressions 
of these counting functions obtained in Section 4.3, we just explain the main origin of all 
possible singularities, and postpone the complete proofs to [12]. 
Consider first Q(0, 0, z). 

Proposition 4.4. The smallest positive singularity of Q(0,0, z) is 

z g = mf{z > : y 2 (z) = y 3 (z)}. (4.6) 

Remark 4.5. We choose to denote the singularity above by Zg ; CIS 3-11 alternative definition 
of z g could be the following: the smallest positive value of z for which the genus of the 
algebraic curve {(x,y) £ C 2 : K(x,y,z) = 0} jumps from 1 to 0. See also Proposition 4.6 
for other characterizations of z g . 

Sketch of the proof of Proposition 4-4- To find the singularities of Q(0, 0, z), we start from 
the expressions obtained in Section 4.3, especially (4.2) and (4.3). 

Consider first the case of the group of order 4. Then the quartic curve Ai z is a circle for 
any z (see Section 4.1), and w(t,z) is a rational function of valuation 1 (see Section 4.2). 
In fact, the singularities come from Yo(t,z), since the branch points xg{z) appear in the 
expression of this function, see (2.4) and (2.6). Hence, the first singularity of Q(0, 0, z) is 
exactly the smallest singularity of the xg(z), and this corresponds to equation (4.6), see 
Proposition 4.6 below. 

Consider now all the remaining cases (i.e., an infinite group, or a finite group of order 
strictly larger than 4). Then the singularity z g appears not only for the same reasons as 
above, but also on account of the CGF w(t,z). Indeed, for z S (0, z g ), the curve M. z 
is smooth, while for z = z g one can show that A4 Z has a non-smooth double point at 
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X(y2(z), z), see [9, 11]. Accordingly, w(t, z) has a singularity at z g , and the behavior of 
w(t, z) in the neighborhood of z g is strongly related to the angle between the two tangents 
at the double point of the quartic curve, see [11]. ■ 

In Proposition 4.4 and in Remark 4.5, we gave two different characterizations of z g . We 
present hereafter five other ones. 

Proposition 4.6. The value of the first singularity z g , introduced in (4.6), can also be 
characterized by the five following equivalent statements. 

(i) z g = inf{z > : x 2 (z) = 2:3(2)}. 

(ii) z g is the smallest positive singularity of the branch points xg(z). 
(Hi) z g is the smallest positive singularity of the branch points yi{z). 

(iv) Zg is the smallest positive root of the resultant (taken w.r.t. x) of the polynomial 
d(x, z). 

(v) z g is related to the minimizer of the Laplace transform of the Sij 's as follows. Define 
(a, (3) as the unique solution in (0, oo) 2 of 

''V 1 ' = °> jSija'P = 0- (4-7) 

Then 

Zg = y 1 s- .<*& ■ m 

Before proving Proposition 4.6, we notice at once that z g (fortunately!) coincides with the 
smallest positive singularity of Q(0, 0, z) found in [7, Section 1.5], which does match (4.7) 
and (4.8). 

Secondly, an easy consequence of any of the five points of Proposition 4.6 is that z g is 
algebraic. Indeed, point (iv) implies that this degree of algebraicity is at most 7, and can 
in fact be strictly smaller than 7, as implied by the following corollary of Proposition 4.6 
(v). 

Corollary 4.7. We have z g = 1/|<S| if and only if Y2-i<i<i = S-i<jXiJ^,i = (this 
happens for 14 models, according to the classification in [5]). Otherwise, z g > 1/|<S|. 

Proof of Proposition 4-6. Only point (v) has to be shown, since the others are easy by- 
products of the definition (4.6) of z g . For this purpose, we use (see (4.6) and Proposition 
4.6 (i)) the equalities 

z g = inf{z > : y 2 (z) = y^{z)} = inf{z > : x 2 (z) = x 3 (z)}. 

Then we have (see, e.g., [9, Chapter 6] or [15, Section 5]) 

X (y2(Zg), Zg) = X 2 (Zg), Y(x 2 {Zg), Zg) = J/ 2 (%)- 
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Hence, the pair (a, (3) = (x2(z g ), 2/2 (zg)) satisfies the system 

d d 
—K{a,(3,z g ) = —K(a,P,z g ), 

which in turn yields (4.7), while (4.8) is a direct consequence of K(a,/3,z g ) = 0. Finally, 
to see why (4.7) has a unique solution in (0, oo) 2 , we refer for instance to [15, Section 5]. 



In [12], while proving Proposition 4.4, we shall besides be able to obtain the precise 
behavior of Q(0, 0, z) near z g . To this end, we use and extend results of [11, 21] to derive 
some fine properties of the CGF, when the Riemann surface passes from genus 1 to genus 
0. 

As for the singularities of (5(1,0, z) and Q(0, l,z), a key point is to locate the point 1 
with respect to the domains &(Ai z ) and Jz?(£ 2 ). Indeed, the expressions of Q(l, 0, z) and 
Q(0, 1, z) written in Section 4.3 depends on this location. 

If 1 belongs to the sets above for any z G (0, z g ), then Q(l, 0, z) and Q(0, 1, z) are simply 
obtained by evaluating the integrals in Proposition 4.3 at x = 1 and y = 1. The first 
singularity of these functions is then z g again. 

Assume now that for some z £ (0, z g ), 1 does not belong to ^(M z ). Then we use (4.4) to 
define Q(1,0, z). Then the singularities of Q(l,0,z) have to be sought among z g and the 
singularities of Yq(1,z) and Xq(Yq(1, z), z). But Xq(Yq(1, z), z) is either equal to 1, or to 
c(Yq(1, z))/a(Yo(l, z)), see [9, Corollary 5.3.5], and accordingly is either regular or has the 
same singularities as Yq(1,z). 

Equations (2.3), (2.4) and (2.6) imply that the singularities of Yq(1,z) necessarily satisfy 
d(l, z) = 0. As a consequence, the smallest positive singularity of Yq(1, z) is given by 

1 

ZY ~ 6(1) + 27^(1)' 

Lemma 4.8. zy e [l/|5|,z 9 ]. 



Proof. The inequality zy ^ 1/|5| is a consequence of the relations 



2y/a(l)c(l) < a(l) + c(l), a(l) + b(l) + c(l) = \S\. 

Also, zy satisfies d(l,zy) = 0. In other words, zy is the smallest positive value of z such 
that 1 is a root of d(l,z). But we have £2(0) = 0, ^3(0) = 00, and X2(z g ) = xs(z g ), see 
Proposition 4.6. Hence, by a continuity argument, one of the points X2(z) or x^(z) must 
become 1 before reaching the other branch point. ■ 

Similarly, we obtain that Q(0, l,z) has a singularity at z g , and possibly at 

1 

zx 



6(1) + 2^5(1)0(1) 
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Concerning the algebraicity of zy and zx, we simply note that these numbers are either 
rational or algebraic of degree two, if and only if a(l)c(l) and a(l)c(l) are not square 
numbers. 

As for Q(l, 1, z), once we know the singularities of Q(0, 0, z), Q(l,Q,z) and Q(0, l,z), it 
is immediate to compute those of Q(l, 1, z) by means of equation (4.5). 

Remark 4.9 (Conclusion). As we have seen, a consequence of the various results of Section 
4 is that the smallest positive singularities of Q(0, 0, z), Q(l, 0, z), Q(0, 1, z) and Q(l, 1, z) 
are algebraic, and sometimes even rational, for all models of walks. 

Acknowledgments. This paper is dedicated to the memory of Philippe Flajolet, who was 
deeply interested in these problems, about which we had together enjoyable preliminary 
discussions. This brutal stroke of fate prevented our friend Philippe of becoming the third 
co-author. . . 
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